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The strong Coulomb interaction between massless Dirac fermions can drive a semimetal-insulator transition in single-layer graphene 
by dynamically generating an excitonic fermion gap. There is a critical interaction strength A c that separates the semimetal phase 
from the insulator phase. We calculate the specific heat and susceptibility of the system and show that they exhibit distinct behaviors 
in the semimetal and insulator phases. 
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1. Introduction 

The low-energy properties of graphene have been widely in- 
vestigated theoretically and experimentally in recent years |Q]]. 
It is well-known that the low-lying elementary excitations of 
graphene are massless Dirac fermions, which have linear dis- 
persion and display quite different behaviors from ordinary elec- 
trons with parabolic dispersion. At half-filling state, the density 
of states of massless Dirac fermions vanishes linearly with en- 
ergy u> near the Fermi level. Due to this fact, there is essen- 
tially no screening on the Coulomb interaction between Dirac 
fermions. The unscreened, long-range Coulomb interaction was 
argued I2 UI0I1 to be responsible for a plenty of unusual physical 
properties, including the logarithmic velocity renormalization 
0,1a], the logarithmic specific heat correction |01, the presence 
of quantum critical point jjfl lg], and the marginal Fermi liquid 
quasiparticle lifetime |9l [lOO . 

When the unscreened Coulomb interaction is sufficiently 
strong, the semimetal ground state of graphene may no longer 
be stable. There exists an interesting possibility that the mass- 
less Dirac quasiparticles and quasiholes are bound into pairs 
through the attractive Coulomb interaction between them. As 
a consequence, the massless Dirac fermions acquire a finite 
mass and the ground state of graphene becomes insulating. This 
semimetal-insulator transition is usually called excitonic insta- 
bility in the literature yl |4j, Illl4l4ll . It can be identified as 



the non-perturbative phenomenon of dynamical chiral symme- 
try breaking conventionally studied in the context of particle 
physics B15U17I1 . Both Dyson-Schwinger equation OLkl ll ill 1411 
and lattice simulation approaches 11121 1 1311 found that such ex- 
citonic instability occurs only when the fermion flavor is less 
than a critical value N c and the Coulomb strength is larger than 
a critical value A c . If we fix the physical fermion flavor N — 2, 
then the semimetal-insulator transition happens at a single crit- 
ical point A c . 

The effective coupling parameter of Coulomb interaction 
can be defined as A = Ng 2 /16 = Ne 2 /(16eohvf ) with eo being 



the dielectric constant and Vp being the effective velocity. In the 
clean limit, the physical magnitude of this parameter is around 
3 or 4 for graphene in vacuum. In the same limit, we found by 
solving gap equation to the leading order of 1 /N expansion that 
the critical strength A c a 1.85 |l4jl . Using analogous gap equa- 
tion approach, the critical coupling is found to be A c « 2.08 and 
A c w 4.16 respectively in Ref. Jj and Ref. fl. In addition, 
the Monte Carlo study lfl3ll performed in lattice field theories 
found that the critical strength A c * 1.74 at N - 2. Our critical 
coupling is much more closer in magnitude to that of Monte 
Carlo study. 

Once a fermion mass gap is generated, the low-energy prop- 
erties of graphene fundamentally change. Below the energy 
scale set by the fermion gap, the density of states of fermions is 
substantially suppressed, which would produce important con- 
sequences. It is interesting to study some observable physi- 
cal quantities those can serve as signatures for the existence 
of excitonic instability. The effects of dynamical fermion gap 
have been discussed by several authors lfl8l - l2011 . In this pa- 
per, we calculate the specific heat and susceptibility of Dirac 
fermions and other low-energy excitations in both semimetal 
and excitonic insulator phases. These quantities can be com- 
pared with experimental results and hence may help to under- 
stand the physical consequence of excitonic instability. 

In the semimetal phase, the Coulomb interaction is not strong 
enough to trigger excitonic pairing instability, but it is strong 
enough to produce unusual properties. As found by Vafek 10], 
the long-range Coulomb interaction gives rise to logarithmic 
T-dependence of fermion specific heat, which is clearly not be- 
havior of normal Fermi liquid. In this paper, we re-derive the 
same qualitative result by a different method. We also calcu- 
late the susceptibility of massless Dirac fermions and show that 
it also exhibits logarithmic T-dependence due to long-range 
Coulomb interaction. 

In the insulator phase, the fermion density of states is sup- 
pressed by the excitonic gap. Intuitively, the specific heat and 
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susceptibility of Dirac fermions should drop significantly from 
their corresponding magnitudes in the semimetal phase. Our 
explicit computations will show that this is true. However, the 
massive Dirac fermions are not the true low-lying elementary 
excitations in the insulating state. At the low energy regime, 
the only degree of freedom is the massless Goldstone boson 
which originates from the dynamical breaking of continuous 
chiral symmetry. The Goldstone bosons make dominant contri- 
bution to the total specific heat at low temperature, but make no 
contribution to the total susceptibility. 

In section [2] we set up the Hamiltonian of the system and 
define the physical quantities in which we are interested. We 
calculate the free energy, specific heat and susceptibility in sec- 
tion[3] We finally summarize the results and discuss some rele- 
vant problems in section|4] 

2. Model and Definitions 

The Hamiltonian of massless Dirac fermions in single layer 
graphene is given by 



Ho 



He = 



N 

<r=l 
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xi/*v(f',r')y iMf',r'), 
where the Coulomb interaction potential J3l is 
, ' da> d 2 k 

U c (t,r) 
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(2) 



(3) 



2tt 2n |k| + U(co, k) ' 

where g 2 = e 2 /eohvf. As mentioned in Introduction, it is con- 
venient to define a dimensionless Coulomb coupling as A = 
Ng 2 116. Usually, Dirac fermion in two spatial dimensions is 
described by two-component spinor field whose 2x2 represen- 
tation can be formulated by Pauli matrices y M = (cr 2 , /<x 3 , icr\). 
However, it is not possible to define a 2 x 2 matrix that anticom- 
mutes with all these matrices. Therefore, there is no chiral sym- 
metry in this representation. Here, we adopt four-component 
spinor field if/ to describe the massless Dirac fermion llj, 16ll . 
The conjugate spinor field is defined as ij/ = ip^yo- The 4x4 
y-matrices can be defined as y M = (0-3, icr\, icri) ® o"3, which 
satisfy the standard Clifford algebra {y^, y v ) - 2g flv with metric 
gftv - diag(l, -1,-1). Obviously, there are two 4x4 matrices 
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which anticommute with all y^. The total Hamiltonian pre- 
serves a continuous U(2N) chiral symmetry tf/ — > e' ani if/. The 
mass term generated by excitonic pairing will break this global 
chiral symmetry dynamically to subgroup U (N) x U (N). Mean- 
while, according to the Goldstone theorem, there appear mass- 
less Goldstone bosons due to the breaking of continuous chiral 
symmetry. These bosons are the only gapless excitations in the 
symmetry broken phase and hence play an important role in de- 
termining the low-energy behaviors of the system. Although 



the physical fermion flavor is actually N - 2, in the following 
we consider a general N in order to perform 1 /N expansion. 
For convenience, we work in units where h = ks = Vf = 1 
throughout the paper. 

The electronic structure of graphene is very special in that 
the ^-conduction bands and 7r*-valence bands touch at two in- 
equivalent K points. This is the reason why the low-energy 
fermionic excitations have a linear dispersion. When the strong, 
long-range Coulomb interaction opens an excitonic gap at the 
Dirac point, the chiral symmetry of total Hamiltonian is broken, 
resembling the non-perturbative phenomenon of dynamical chi- 
ral symmetry breaking in QED3 1 17]. This mechanism was first 
proposed in graphene by Khveshchenko yl and has been ex- 
tensively studied flE -14] in the following years. 

Gusynin et al. discussed the influence of excitonic fermion 
gap on various transport quantities, including electrical and Hall 
conductivity (H, 18]. The results were compared directly with 
the experiments in graphene. Recently, Kotov et al. studied the 
effect of fermion gap on the interacting potential [ 19] and found 
an effective weak confinement of fermions. They also argued 
that the massive phase exhibits much more interesting behav- 
ior than the massless one. The effect of fermion gap on quasi- 
particle lifetime and spectral function was discussed in lEoll . 
This kind of excitonic instability may also exist in other cor- 
related electron systems than graphene. For instance, it was 
suggested by one of the authors that such instability can pro- 
vide a qualitative understanding on the field-induced thermal 
metal-insulator transition observed in the vortex state of high 
temperature cuprate superconductor 112 111 . 

In this paper, we calculate the specific heat and susceptibil- 
ity by including the effect of Coulomb interaction in both the 
semimetal and insulator phases. These are physical quantities 
those can be measured by experiments and hence can help us 
to build interesting connections between theoretical predictions 
and experimental facts. Technically, we will follow the proce- 
dures utilized in the paper of Kaul and Sachdev l22ll . In this 
framework, all propagators and correlation functions are writ- 
ten in the Matsubara imaginary time formalism. 

At finite temperature, the fermion propagator is 



G(ia>„, k) 



1 



iu>„y -y-k-m 



(4) 



where a>„ = (2n + 1)7tT is the fermion frequency. Although 
generally the excitonic fermion gap should depend on momen- 
tum, energy, and temperature, we assume a constant mass gap 
m throughout the paper to simplify calculations. 
The bare Coulomb interaction function is simply 



g 2 * 
D (co m ,q) = 2- = — . 

" n 



2q 



(5) 



In an interacting electron gas, the collective excitations screen 
the bare Coulomb interaction and convert Do to 



D(w,„,q, T) 



1 
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f a + «n( Wm ,q,r) 
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where the polarization function n is defined as 

' 2 k Tr[ r ofro« + *)] 



U(cj m ,q, T) = -NT 
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(7) 



with <7o = <Jj m — 2mnT and Icq = oj n — (2n + \)nT . 

The first two orders in 1 /N expansion of free energy T are 

r = Nf f +f f , (8) 

where f°f is the leading, noninteracting term and f 1 ^ the cor- 
rections from the Coulomb interaction. The leading term of 
fermion free energy is defined as ff = T%f ^ m[G(iW„,k)]. 

d), } 

To calculate the free energy, we will first sum over the Matsub- 
ara frequencies w„ and then perform the integration over the in- 
termediate variables and momentum k, dropping all terms those 
are independent of temperature and volume B23I1 . The volume 
factor is neglected throughout this paper and we only consider 
free energy in unit volume. The interaction correction to free 
energy is given by 



T(T) = T 



v f£i 

(2* 



q yiMD- 1 ]. 



(9) 



We choose the zero temperature free energy T(T = 0) as the 
reference free energy 10, Hi)], and then define the following reg- 
ularized free energy 

ff = T{T) 



0) 
In 



D-\aj m ,q, T) 



D-Hu> m ,q,T = 0) 



(10) 



Here, we follow the strategy of Ref. 1I22I1 and introduce a mag- 
netic field H. For fermions, the field shifts frequency as a>„ — > 
co„ - OH, where 8 — ± 1 . The specific heat Cy and susceptibility 
Xf can be defined as 

„d 2 T 



C 



Xf 



dT 2 



JVC?/ + C 
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(12) 



dH 2 \h=o 

which are divided to free and interaction terms, respectively. 
For a normal Fermi liquid, the specific heat and susceptibility 
should behave as Cy °c T 2 andxf 00 T according to the analysis 
in Ref. H22l 12511 . If we write the specific heat as Cy = ^lc v T 2 , 
then Mc v should be 



Similarly, the susceptibility can also be written as Xf 
with Zh x , being 



'Xs 



if 



(13) 



(14) 



In the presence of fermion mass m, the specific heat (suscepti- 
bility) no longer behaves as oc T 2 (oc T). However, in order to 
make direct comparison, we still express specific heat (suscep- 
tibility) in terms of Mc v i&xi)' which w iH depend on tempera- 
ture T. The definitions presented in this section will be used to 
calculate the free energy, specific heat, and susceptibility in the 
next section. 



3. Specific heat and susceptibility 

3.1. Leading terms 

In the presence of a constant fermion mass m, the noninter- 
acting free energy is 

/d 2 \i 
— 5 ln[G(kfc,k)] 
(2n) 2 

= -T f ^ 

J 2n 2 
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(15) 



Here, Li2 and Li3 are polylogarithmic functions. It is easy to 
get the following noninteracting term for fermion specific heat 
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(16) 



This function is plotted in Fig. |l(a)| Taking the m = limit of 
Eq. (TToT ). the specific heat in the semimetal phase is 



r o/ _ 9£(3) 2 
v ~ 2n ' 



(17) 



with M J = Similarly, the noninteracting susceptibility is 
given as 



(18) 



which is plotted in Fig. |l(b)| Taking the m = limit of Eq. 
( fT8l . the susceptibility in the semimetal phase is 



Xf = T, 

1 n 



(19) 



withjfl^ = 

From Eq. ( [ToT l. Eq. (IT~8b . and Fig. Q~l it is easy to see that the 
T-dependence of specific heat and susceptibility in the insulat- 
ing phase differs significantly from the corresponding oc T 2 and 
oc T behaviors in the semimetal phase. This is not unexpected 
because the excitonic gap strongly suppresses the fermionic ex- 
citations at low temperature. 

However, although the low-energy fermion excitations are 
strongly suppressed in the insulating phase, there exists another 
kind of gapless excitation: Goldstone boson. The presence of 
gapless Goldstone bosons is the characteristic property of exci- 
tonic instability. They are composed of Dirac fermions (quasi- 
particles) and anti-fermions (quasiholes), but carry no electric 
charge themselves. The Goldstone bosons do not contribute to 
the susceptibility because they do not couple to external mag- 
netic field H, but they do contribute to the total specific heat of 
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the system. In particular, the free energy of Goldstone bosons 
is 



2^(3)^ 



while the corresponding specific heat is 



12^7*. 



(20) 



(21) 



which has been obtained in Ref . II26I1 . 

Apparently, a oc T 2 term of specific heat appears in both 
the semimetal phase and the insulating phase. It seems that 
these two phases have similar specific heat, albeit contributed 
from different elementary excitations. However, such similarity 
actually does not exist because it disappears once the interaction 
correction to free energy is incorporated. 
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Figure 1 : (a) The dependence of fermion specific heat on fermion gap j . The 



a"./ 



%0) 



at m = 0. 



top dashed line corresponds to the non-interacting term - -jj- 
(b) The dependence of fermion susceptibility on fermion gap j . The top dashed 

line corresponds to the non-interacting term S$£ = ^= at m = 0. It is obvious 
that fermion gap suppresses fermion specific heat and susceptibility strongly. 



3.2. Interaction corrections 

We now include the interaction correction to the free en- 
ergy. Note the Goldstone bosons are neutral, so the Coulomb 
interaction only affects the free energy of Dirac fermions. To 
calculate the free energy f x * , we should first know the polar- 
ization function. In the presence of finite fermion mass m and 



external magnetic field H, the polarization function YlUom, q, T) 
can be calculated by the methods presented in fl4lr26l, 271. Here 
we only write down the final expression: 

U(a> m , q, T) 



! (1 - 2x)(jj m sin(ixw,„ + ^6H) 



nt r 

= dxln[4D m ] 

n Jo 

I dx 

4tt Jo 

X 



N f 1 , m 2 + x(l - x)col sinhifEf) 

— I dx , 

2/r Jo E f D m 



(22) 



at finite temperature. Here we introduced the following abbre- 
viated notations 



Ef = ^m 2 + x{\ - x){q 2 + oj 2 m ), 



,2 Ef, 



, xco m + 6H 



D m = cosh z (— ) -sin z (- 



-)• 



(23) 
(24) 



The zero-temperature limit of polarization function (Eq. (l22l ) 



is 



N f 1 x(l-x)q 
n(w m ,q,T = 0) = - dx ;H 
x Jo E f 



(25) 



Using Eq. (l22l and Eq. ( 1251 ). the free energy f 1 ? (Eq. (flOl l) can 
be directly computed. 

In the semimetal phase with m = 0, the polarization func- 
tion is 



n(ai m , q, T,m-0) 



q 2 NT 



q 2 + < 



it r 

n Jo 



dx 



xln 



2 1 cosh(-^) + cos(— — + 0—) 



at finite temperature with £/o = s/xjl - x)(q 2 + w 2 ) and 



n(w,„,q,r = 0,m = 0) = 



N 



8 Vq 2 + <4 ' 



(26) 



(27) 



at zero temperature. Using these expressions, the free energy of 
Dirac fermion is written as 

4r 3 r°° r 1 

f lJ = ^ J g J o d y J o dx{K(x,y,q) + J(x,y,q)},m 
«p[-2qVx(l-x)(l+y a )] 



where 
K(x,y,q) 



J(x,y, q) 



tt(1 +y 2 ) 


1+ r^} 




2exp[-qV^(l -x)(l +y 2 ) 


n(l +y 2 ) 






x cos(xyq + 1 







(29) 



(30) 



Here, a variable y = ^ is introduced, with a> being the contin- 
uous form of u> m when T — > 0. For finite y, K(x, y, q) damps 
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rapidly with growing y, so y ~ makes the dominant contribu- 
tion to the free energy. We can expand the function K(x, y, q) 
near this point and obtain 



fl f = 



4-T 3 C7 r°° r 1 

— 5- q^q dy dxK(x,y,q) 
x Js Jo Jo 

— 5- I qifq dy I dx 

x Js Jo Jo 

exp[-qV(l-4x 2 )(l+y 2 )] 



(1 + y 2 ) 



8T 3 „_ 

— qdq 

" Js Jo 



cfy- 



(1+y 2 ) 



r 



exp[2qVlV0^7)] 



where 



tan 



riW = 1 + 



7T 

2I' 



(31) 



(32) 



Comparing with K(x, y, q), the form of J(x, y, q) is more com- 
plicated owing to the cosine term cos{xyq+8H/T). The compu- 
tation becomes difficult if we make Taylor expansion of the co- 
sine function. By plotting the dependence of function J(x, y, q) 
on its variables, we found that the dominant regime is x ~ 
0, y ~ 0. Hence, we simply take J(x, y, q) as 



J(x,y,q) 



2exp[-q +y 2 )] H 



(1+y 2 ) 









cos(0~), 



(33) 



which then leads to 

-3 



4 T i rr r°= r l 
f 2 = — J Q^q J dy J dxJ(x,y,q) 

T T 

Taking H — 0, the total free energy now has the form 



(34) 



f lf = —r!(A)T 3 ln^. 



20 



A 



T 

It is easy to get the following specific heat and susceptibility 



(35) 



1/ 



C 



1/ 



120 , A 

7T-' / 



(36) 
(37) 



Here, the ultraviolet cutoff A can be taken to be of order lOeV, 
which is determined by ~ cr 1 with lattice constant a = 2.46A. 



From these results, we know that both specific heat and sus- 
ceptibility of massless Dirac fermions exhibit logarithmic In- 
dependence due to long-range Coulomb interaction. These are 
non-Fermi liquid behaviors. 

The appearance of such singular fermion specific heat was 
first pointed out by Vafek fj]. Here, we obtained the same qual- 
itative T-dependence by a different method. In Ref. [7], the 
calculation of free energy was performed on the basis of the 
retarded vacuum polarization functions and retarded fermion 
propagator G Kt (u>, k) = 1 k , while in our case the polar- 
ization functions and fermion propagator are expressed in the 
Matsubara formalism. Strictly speaking, these two polariza- 
tion functions are equivalent and should lead to the same re- 
sults. We numerically compute the free energy using both the 
polarization functions obtained in the present paper and that 
in Ref. |7J], and found that the results are very close to each 
other (the maximum proportional error of the coefficient ^ x 
is < 5%). In order to get an analytic expression for free en- 
ergy, some approximations to the polarization functions is un- 
avoidable. In Ref. [7], the dominant contribution of polar- 
ization function comes from y = ^ ~ 1 at both y > 1 and 
y < 1 regions (after analytic continuation the momentum be- 
comes q = ^q 2 - {a> + id) 2 = q yjl — y 2 ), while in our calcula- 



tion the dominant momentum region isy = =* ~ (y = — 







in the continuous form). For this reason, our analytic expres- 
sion for the free energy differs from that of Ref. |0] (the ap- 
proximation of J(x,y,q) might partly explain the difference). 
After comparing the analytical results with numerical results, 
we found that our analytical result is slightly lower than the nu- 
merical result while the analytical result in Ref. [7] is slightly 
greater than the numerical result. For A = 1, the analytical and 
numerical results for the coefficient are 0.277 and 0.23 1 respec- 
tively in our work and 0.200 and 0.225 respectively in Ref. 

The contribution of Coulomb interaction to the free energy 
in the semimetal phase with m — is shown in Fig. [2] The 
free energy behaves as oc T 3 In T (logarithmic correction) for 
several different values of A. 

We now turn to the insulator phase where m + 0. The free 
energy can be obtained by substituting Eq. (1221 and Eq. ( 1251 ) 
into Eq. ( [Tot . The dependencies of specific heat and suscepti- 
bility on different fermion mass m for A — 4 are shown in Fig. 
|3(a)| and |3(b)| respectively. The results for other choices of A 
are similar and thus not shown. Here, we use the absolute val- 
ues |^| = -Jl l c f y and |^| = -m\ f r instead of ^ and jC x f f 
which are negative. From Fig. |3(a)| and |3(b)| we see that the 
fermion gap leads to remarkable suppression of the interaction 
correction to fermion specific heat and susceptibility. 

In summary, in the semimetal phase the long-range Coulomb 
interaction gives rise to non-Fermi liquid behavior of specific 
heat and susceptibility. In the insulator phase, the fermion spe- 
cific heat and susceptibility are both significantly suppressed 
by the excitonic gap, but the total specific heat has a finite value 
due to the massless Goldstone bosons. 
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Figure 2: Interaction correction to free energy in the semimetal phase with 
m = for different interaction strength A. The red line is the reference free 
energy ln£. It appears that free energy displays the same logarithmic behavior 
for different A. 



4. Conclusion and discussion 

In this paper, we calculated the specific heat and suscepti- 
bility in graphene. The ground state of graphene is semimetal 
when the Coulomb interaction strength A < A c , but becomes in- 
sulator when A > A c . The most prominent feature of semimetal 
phase is the appearance of logarithmic T-dependence of spe- 
cific heat and susceptibility due to long-range Coulomb inter- 
action. These are non-Fermi liquid behaviors. In the insulating 
phase, because the interaction correction to fermion excitations 
is strongly suppressed by the excitonic gap, the total specific 
heat is solely determined by the contribution from Goldstone 
bosons, while the susceptibility drops significantly. Apparently, 
both specific heat and susceptibility manifest quite different be- 
haviors in the two sides of the critical point A c . 

Note that the semimetal and insulator phases both contain 
massless excitations: massless Dirac fermion in the former and 
massless Goldstone boson in the latter. They have different 
statistics and exhibit completely different behaviors. For ex- 
ample, the massless Dirac fermions can transfer heat current 
and produce a universal thermal conductivity lEii at T = 0, 
while the Goldstone bosons only contribute a oc r 3 term, which 
vanishes rapidly as T — > 0. The massless Dirac fermions also 
gives rise to a universal electric conductivity B28I1 . although the 
predicted electronic conductivity is at invariance with experi- 
mental result (the famous missing n). The Goldstone bosons do 
not contribute to electric conductivity since they are neutral. 

We should point out that the Goldstone bosons are exactly 
massless only when the Lagrangian respects a continuous chiral 
symmetry. If the continuous chiral symmetry is explicitly bro- 
ken by some contact four-fermion interaction, then the Gold- 
stone bosons are no longer strictly massless. Instead, they have 
a small mass as the result of dynamical breaking of appropriate 
continuous chiral symmetry [29]. In this case, our discussion 
and calculation about the free energy contribution from Gold- 
stone bosons should be modified and a small mass should be 
included. In reality, there are various four-fermion interactions 
in the graphene 130113 ill . If the contact four-fermion interaction 
has the form oc (ifryoip) 2 , then the continuous chiral symme- 
try is not explicitly broken and the Goldstone bosons are still 




(b) 



Figure 3: (a) Interaction correction to fermion specific heat in the insulator 
phase with m + at A = 4. (b) Interaction correction to susceptibility in the 
insulator phase with m # at A = 4. The suppressing effect of fermion gap is 
apparent. 



massless. If the four-fermion interaction term is oc (tfri//) 2 , then 
the system has only discrete chiral symmetry and there are no 
massless Goldstone bosons [14, l3~ltl . Therefore, the specific 
heat of Goldstone bosons presented in Sec[3]is valid only when 
the continuous chiral symmetry is not explicitly broken by any 
four-fermion interaction term. 

We finally comment on the validity of 1 /N expansion. The 
excitonic insulating transition requires the Coulomb interaction 
between Dirac fermions be sufficiently strong. In this strong 
coupling regime, 1 /N seems to be the only available expansion 
parameter, even if it is not small (N — 2 for graphene). In our 
specific case, the fermion mass plays the dominant role in the 
insulator phase. It suppresses significantly the Coulomb inter- 
action contribution to fermion specific heat. This implies that, 
within the 1 /N expansion, the next-to-leading order contribu- 
tion could be neglected since it is much less than the leading 
order contribution. It is reasonable to speculate that higher or- 
der corrections in l/N expansion are also suppressed by the 
dynamical fermion mass. In the semimetal phase, there is no 
such suppressing effect, so higher order corrections might be 
important. As shown in the context, the analytical calculation 
of next-to-leading order correction is already very complicated, 
including higher order corrections will make analytical calcu- 
lation intractable. The specific heat of massless Dirac fermions 
may be analyzed by renormalization group approach |0, l32ll . 
which found power-law T@ behavior after summing up all or- 
ders of logarithmic corrections ll6i [32Tl . However, the exponent 
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/? can only be calculated by performing 1 /N expansion. There- 
fore, the validity of 1 /N expansion also needs to be studied in 
this approach. 
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